Transcriptome differential expression analysis of defoliation in different lemon varieties under drought treatment.

'Allen Eureka' is a bud variety of Eureka lemon with excellent fruiting traits, but severe winter defoliation affects the following year's yield, and the response mechanism of lemon defoliation is currently unknown. Two lemon cultivars ('Allen Eureka' and 'Yunning No. 1') with different defoliation traits were used as materials to investigate the molecular regulatory mechanisms of different leaf abscission periods in lemons. The petiole abscission zone was collected at three different defoliation stages, namely, the predefoliation stage (k15), the middefoliation stage (k30), and the postdefoliation stage (k45). Transcriptome sequencing was performed to analyze the gene expression differences between these two cultivars. A total of 1141, 2695, and 1433 differentially expressed genes (DEGs) were obtained in k15, k30, and k45, respectively, and the number of DEGs in k30 was the largest. GO analysis revealed that the DEGs between the two cultivars were mainly enriched in processes related to hydrolase activity, chitinase activity, oxidoreductase activity, and transcription regulator activity in the defoliation stages. KEGG analysis showed that the DEGs were concentrated in k30, which involved plant hormone signal transduction, phenylpropanoid biosynthesis, and biosynthesis of amino acids. The expression trends of some DEGs suggested their roles in regulating defoliation in Lemon. Seven genes were obtained by WGCNA, including sorbitol dehydrogenase (CL9G068822012_alt, CL9G068820012_alt, CL9G068818012_alt), abscisic acid 8'-hydroxylase (CL8G064053012_alt, CL8G064054012_alt), and asparagine synthetase (CL8G065162012_alt, CL8G065151012_alt), suggesting that these genes may be involved in the regulation of lemon leaf abscission.

Transcriptome differential expression analysis of defoliation of two different lemon varieties.

'Allen Eureka' is a bud variety of Eureka lemon with excellent fruiting traits. However, it suffers from severe winter defoliation that leads to a large loss of organic nutrients and seriously affects the tree's growth and development as well as the yield of the following year, and the mechanism of its response to defoliation is still unclear. In order to investigate the molecular regulatory mechanisms of different leaf abscission periods in lemon, two lemon cultivars ('Allen Eureka' and 'Yunning No. 1') with different defoliation traits were used as materials. The petiole abscission zone (AZ) was collected at three different defoliation stages, namely, the pre-defoliation stage (CQ), the mid-defoliation stage (CZ), and the post-defoliation stage (CH). Transcriptome sequencing was performed to analyze the gene expression differences between these two cultivars. A total of 898, 4,856, and 3,126 differentially expressed genes (DEGs) were obtained in CQ, CZ, and CH, respectively, and the number of DEGs in CZ was the largest. GO analysis revealed that the DEGs between the two cultivars were mainly enriched in processes related to oxidoreductase, hydrolase, DNA binding transcription factor, and transcription regulator activity in the defoliation stages. KEGG analysis showed that the DEGs were concentrated in CZ and involved plant hormone signal transduction, phenylpropanoid biosynthesis, glutathione metabolism, and alpha-linolenic acid metabolism. The expression trends of some DEGs suggested their roles in regulating defoliation in lemon. Eight gene families were obtained by combining DEG clustering analysis and weighted gene co-expression network analysis (WGCNA), including ß-glucosidase, AUX/IAA, SAUR, GH3, POD, and WRKY, suggesting that these genes may be involved in the regulation of lemon leaf abscission. The above conclusions enrich the research related to lemon leaf abscission and provide reliable data for the screening of lemon defoliation candidate genes and analysis of defoliation pathways.

Comparative physiochemical and transcriptomic analysis reveals the influences of cross-pollination on ovary and fruit development in pummelo (Citrus maxima).

'Shuijingmiyou' pummelo (SJ), one of the most popular fruits in Yunnan province of China, is of relatively low fruit shape (FS) quality. In this study, we compared the FS promoting effects of cross pollinations using pollens from seven pummelo varieties, and found that 'Guanximiyou' pummelo (GX) cross-pollination showed the best FS promoting effects on SJ fruits by shortening its fruit neck. To explore the underlying mechanism, physiochemical and transcriptomic differences between self- and cross-pollinated SJ ovaries (SJO and GXO) were investigated. Higher salicylic acid, gibberellin and indole acetic acid contents and superoxide dismutase, peroxidase and catalase activities, and lower polyphenol oxidase activity were determined in GXO compared with SJO. Enrichment analysis of the identified 578 differentially expressed genes (123 up-regulated and 455 down-regulated) in GXO showed that genes involved in solute transport, RNA biosynthesis, phytohormone action and cell wall organization were significantly enriched. The results obtained in this study will be helpful in understanding the influences of cross-pollination on pummelo ovary and fruit development, and can provide the basis for clarifying the underlying mechanism of cross-pollination improved fruit quality.

Pangenome analysis provides insight into the evolution of the orange subfamily and a key gene for citric acid accumulation in citrus fruits.

The orange subfamily (Aurantioideae) contains several Citrus species cultivated worldwide, such as sweet orange and lemon. The origin of Citrus species has long been debated and less is known about the Aurantioideae. Here, we compiled the genome sequences of 314 accessions, de novo assembled the genomes of 12 species and constructed a graph-based pangenome for Aurantioideae. Our analysis indicates that the ancient Indian Plate is the ancestral area for Citrus-related genera and that South Central China is the primary center of origin of the Citrus genus. We found substantial variations in the sequence and expression of the PH4 gene in Citrus relative to Citrus-related genera. Gene editing and biochemical experiments demonstrate a central role for PH4 in the accumulation of citric acid in citrus fruits. This study provides insights into the origin and evolution of the orange subfamily and a regulatory mechanism underpinning the evolution of fruit taste.

[Fruit variation and geographical distribution of citron].

The ripe dried fruit of citron(Citrus medica) is one of the important sources of Chinese herb Citri Fructus. At the same time, it is also grown for edible and ornamental uses. There are many species and abundant genetic variation. To clarify the intraspecific variation and resource distribution of citron, this study investigated the variation in 11 citron fruits, basically covering the main species in China, including Xiaoguo citron(C. medica var. ethrog), Goucheng(C. medica var. yunnanensis), Muli citron(C.medica var. muliensis), Dehong citron(C.medica×Citrus spp.), Fuzhou citron(C.medica×C.grandis?), Mawu(C.medica×C.grandis?), Cangyuan citron, Binchuan citron, Sweet citron, Big citron, and Small citron. The natural communities of citron were proved to be mainly distributed in the southwestern and western Yunnan and southeastern Tibet of China, with Yunnan, Sichuan, Guangxi, Chongqing, Hubei, and Zhejiang identified as the main production areas. Citron has also been widely grown in India, the Mediterranean region, and the Caribbean coast countries. The field investigation revealed the large-scale intraspecific variation of citron fruits. Most of the fruits are oval-like or sphere-like in shape. The fruits are green when raw and yellow when ripe, with oil cell dots on the skin, stripe-likes running from top to bottom, and bulge at the top. Usually, in the smaller citron fruits, the pulp and juice vesicles are better developed and the central columella is tighter. By contrast, the juice vesicles and central columella in larger fruits became more vacant, with carpels visible, and the apex segregation and development of the carpels is one of the reasons for variation. These variations should be given top priority in the future variety selection and breeding, and the quality differences of different citron species and their mechanisms should be further studied. In particular, variety selection and classification management according to their medicinal or edible purposes will provide scientific and technological supports for the orderly, safe, and effective production of citron products consumed as food and medicine.

Transcriptome and Metabolome Comparison of Smooth and Rough Citrus limon L. Peels Grown on Same Trees and Harvested in Different Seasons.

Background: Farmers harvest two batches fruits of Lemons (Citrus limon L. Burm. f.) i.e., spring flowering fruit and autumn flowering fruit in dry-hot valley in Yunnan, China. Regular lemons harvested in autumn have smooth skin. However, lemons harvested in spring have rough skin, which makes them less attractive to customers. Furthermore, the rough skin causes a reduction in commodity value and economical losses to farmers. This is a preliminary study that investigates the key transcriptomic and metabolomic differences in peels of lemon fruits (variety Yuning no. 1) harvested 30, 60, 90, 120, and 150 days after flowering from the same trees in different seasons. Results: We identified 5,792, 4,001, 3,148, and 5,287 differentially expressed genes (DEGs) between smooth peel (C) and rough peel (D) 60, 90, 120, and 150 days after flowering, respectively. A total of 1,193 metabolites differentially accumulated (DAM) between D and C. The DEGs and DAMs were enriched in the mitogen-activated protein kinase (MAPK) and plant hormone signaling, terpenoid biosynthesis, flavonoid, and phenylalanine biosynthesis, and ribosome pathways. Predominantly, in the early stages, phytohormonal regulation and signaling were the main driving force for changes in peel surface. Changes in the expression of genes associated with asymmetric cell division were also an important observation. The biosynthesis of terpenoids was possibly reduced in rough peels, while the exclusive expression of cell wall synthesis-related genes could be a possible reason for the thick peel of the rough-skinned lemons. Additionally, cell division, cell number, hypocotyl growth, accumulation of fatty acids, lignans and coumarins- related gene expression, and metabolite accumulation changes were major observations. Conclusion: The rough peels fruit (autumn flowering fruit) and smooth peels fruit (spring flowering fruit) matured on the same trees are possibly due to the differential regulation of asymmetric cell division, cell number regulation, and randomization of hypocotyl growth related genes and the accumulation of terpenoids, flavonoids, fatty acids, lignans, and coumarins. The preliminary results of this study are important for increasing the understanding of peel roughness in lemon and other citrus species.

Characterization of the complete chloroplast genome of 'Yunning No.1' lemon (Citrus limon).

'Yunning No.1' lemon, a mutant of Eureka lemon, is originally found in Yunnan province of China and is the main cultivated lemon variety there. In this study, we assembled and annotated its chloroplast genome using Illumina Hiseq-2500 whole genome re-sequencing data. Its chloroplast genome is 160,141 bp in size, containing a 87,754 bp large single copy region, a 18,385 bp small single copy region and a pair of 27,001 bp inverted repeat region. Like many citrus species, 114 unique genes (including 80 protein-coding genes, 30 tRNAs and 4 rRNAs) could be identified from the chloroplast genome of 'Yunning No.1'. Phylogenetic analysis revealed that the 'Yunning No.1' chloroplast genome was closest to Citrus maxima.

Variations of Flavonoid Composition and Antioxidant Properties among Different Cultivars, Fruit Tissues and Developmental Stages of Citrus Fruits.

A large number of biologically active compounds are present in ripe citrus fruits. However, few studies have been focused on the changes in flavonoids and the evolution of antioxidant activity during citrus fruit growth. In this study, fruits of five citrus cultivars cultivated in China were sampled at 60-210 days post-anthesis (DPA) at intervals of 30 days. The amounts of main flavonoids in the peel and pulp were analyzed by HPLC and their activities were studied by DPPH, ABTS and FRAP. The results showed that the contents of hesperidin, diosmin, eriodictyol, rutin and nobiletin increased before 90 DPA and then decreased with the growth and development of fruits, but an opposite tendency was observed for naringin and narirutin. The antioxidant activities in citrus peel and pulp were found to be significantly correlated with some flavonoids. The results may be of guiding values in citrus production and utilization of citrus fruit by-products.

Characterization of the complete chloroplast genome of Yuanjiang wild Ichang papeda (Citrus cavaleriei) in China.

Ichang papeda (Citrus cavaleriei) is an endemic perennial plant in China. In this study, we assembled and annotated the complete chloroplast genome of Yuanjiang wild Ichang papeda using Illumina Hiseq-2500 sequencing data. The chloroplast genome is constituted of 160,996 bp, containing a 87,634 bp large single-copy region, a 18,762 bp small single-copy region, and a pair of 27,300 bp inverted repeat regions. The chloroplast genome contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and 4 rRNAs. Phylogenetic analysis showed that the relationship between the chloroplast gennomes of C. cavaleriei and C. reticulata is the closest, which consistently support their chloroplast relationships.

Characterization of the complete chloroplast genome of Citrus hongheensis, a key protected wild plant in Yunnan province of China.

Citrus hongheensis is a key protected wild plant endemic to the Honghe river region in southeastern Yunnan, China. In the present study, its chloroplast genome was successfully assembled and annotated based on the Illumina Hiseq-2500 whole genome re-sequencing data. The chloroplast genome is 160,275 bp in size. Its large single copy region, small single copy region and inverted repeat region is 87,886 bp, 18,387 bp and 27,001 bp, respectively. Totally, 114 unique genes, including 80 protein-coding genes, 30 tRNAs and 4 rRNAs, were identified from the C. hongheensis chloroplast genome. According to the phylogenetic analysis result, the relationship between the chloroplast genome of C. hongheensis and C. maxima was found to be the closest.